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CKN theory of singularities of 
weak solutions of the Navier-Stokes equations 

Giovanni Gallavotti 
Fisica, Romal and I.N.F.N. 

Abstract: A series of five lectures delivered at the C.I.M.E. 
course on “Mathematical foundations of turbulent viscous 
flows”, 1-6 September, 2003 

§0. Notations. 

The lectures are devoted to a complete exposition of the theory of sin¬ 
gularities of the Navier Stokes equations solution studied by Leray, in a 
simple geometrical setting in which the fluid is enclosed in a container 
Q with periodic boundary conditions and side size L. The theory is due 
to the work of Scheffer, Caffarelli, Kohn, Nirenberg and is called here 
CKN-theory as it is inspired by the work of the last three authors which 
considerably improved the earlier estimates of Scheffer. 

Although the theory of Leray is well known I recall it here getting at the 
same time a chance at establishing a few notations, [Ga02]. 

(1) An underlined letter , e.g. A, A,. .., denotes a 3-dimensional vector 

( i.e . three real or complex numbers) and underlined partial derivative 
symbol 8, 8 ,... denotes the gradient operator (81,82,83). A vector field 

u is a function on fh 

(2) repeated labels convention is used (labels are letters or other) when 
not ambiguous: hence A- B_ or A ■ B means sum over i of AiBi. Therefore 

8-u, if u is a vector field, is the divergence of u, namely 8iUi. 8-8 = A 
is the Laplace operator. 

(3) Multiple derivatives are tensors, so that 88f is the tensor d tJ f. The 
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1/2 (f 2 ) is the space of the square integrable functions on D: the squared 
norm of / G L2 will be ||/||\ d = \f(x)\ 2 sx. 

( 4 ) The Navier-Stokes equation with regularization parameter A is 

( 0 . 1 ) 

where the unkowns are u(x, t),p(x, t) with zero average and, [Ga 02 ], 

(i) u is a divergenceless field, p is a scalar field 

(ii) (u)\ = f^x a(t — y)v±{y)dy and xa is defined in terms of a C°°(Q) 
function x —> x(t) A 0 not vanishing in a small neighborhood of the origin 
and with integral f xGh) dx = l: the function xGk) can be regarded as a 
periodic function on D or as a function on R 3 with value 0 outside U, as 
we shall imagine that D is centered at the origin, to fix the ideas. For 
A > 1 also the function Xa(t) *= A 3 x(Ax) can be regarded as a periodic 
function on D or as a function on l? 3 : it is an “approximate Dirac’s 5 - 
function”. Usually the NS equation contains a volume force too: here we 
set it equal to 0. 1 

(iii) The initial datum is a divergenceless velocity field u° G I/2(U) with 
0 average; no initial datum for p as p is determined from u°. 

( 5 ) A weak solution of the NS equations with initial datum u 0 G I/2(U) is a 
limit on subsequences of A —> 00 of solutions u x ,p x of ( 0 . 1 ). This means 
that the Fourier transforms of u x , and p x exist and have components 
u^(t) = f ci e~' l —u x (x)dx ( k = ^n, n G Z 3 ), and p^{t) which have a 
limit as A -> 00 (on subsequences) for each k ^ 0 and the limit of the 
u fc (t) is absolutely continuous. This is equivalent to the existence of the 
limits of the L2 products (u(t), <^)l 2 and (p(t),V 0 u l° r 1 £ ( 0 , 00] 
and for all test functions f{x)^{x)- There might be several such limits 
(i.e. the limit may depend on the subsequence) and what follows applies 
to any one among them. 

This is a simplicity assumption as the extension of the theory to cases with time inde¬ 
pendent smooth ( e.g. C °°) volume forces would be immediate and just a notational 
nuisance. 


U 


= oAu — 


U- [U) X 


d u — dp , 


d ■ u = 0 , 


u dx = 0 
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The core of the analysis will deal with the regularized equation and the 
properties of its solutions, which are easily shown to be C°° in x G h 
and in t G (0, oo) if the intial datum is u G The limit A —> oo will 

be taken at the end and it is where the theory becomes non conctructive 
because there is need to consider the limit on subsequences. 

Of course the point is to obtain bounds which are uniform in A —> oo and 
the limit A —> oo only intervenes at the end to formulate the results in a 
nice form. 

The theory of Leray is based on the following a priori bounds, see section 
3.2 in [Ga02], on solutions of (0.1) with initial datum u 0 with L 2 square 
norm Eq 

\\u x (t)\\l < E 0l J dT\\du X ( t)\\1 < i E 0 v~ x (0.2) 

satisfied by the solution u x . 

The notes are extracted from reference [Ga02] to which the reader is 
referred for details on the above results and have been made independent 
from [Ga02] modulo the above results (in fact, essentially, only modulo 
the statements on the regularized equation (0.1)). 

The proof is conceptually quite simple and is based on a few (clever) 
a priori Sobolev inequalities: the estimates are discussed in Sect. 1-3, 
which form an introduction. 

Their application to the analysis of (0.1) is in Sect. 4 where the main the¬ 
orems are discussed and the CKN main result is reduced via the inequal¬ 
ities to Scheffer’s theorem. The method is a kind of multiscale analysis 
which allows us to obtain regularity provided a control quantity, identified 
here as the “local Reynolds numbeE on various scales, is small enough. 
Unfortunately it is not (yet) possible to prove that the local Reynolds 
number is small on all small regions (physically this would mean that in 
such regions the flow would be laminar, hence smooth on small scale). 
However the a priori bounds give the information that the local Reynolds 
number must be small near many points in U and, via standard tech¬ 
niques, an estimate of the dimension of the possibly bad points follows. 
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The application to the fractal dimension bound is essentially an “abstract 
reasoning” consequence of the results of Sect. 4 and is in Sect. 5. 

The proofs of the various Sobolev inequalities necessary to obtain the 
key Scheffer’s theorem and of the new ones studied by CKN is decribed 
concisely but in full detail in the series of problems at the end of the 
text: the hints describe quickly the various steps of the proofs (however 
without skipping any detail, to my knowledge). 

§1. Leray’s solutions and energy. 

The theory of space-time singularities will be partly based upon simple 
general kinematic inequalities , which therefore have little to do with the 
Navier-Stokes equation, and partly they will be based on the local en¬ 
ergy conservation which follows as a consequence of the Navier-Stokes 
equations but it is not equivalent to them. 

Energy conservation for the regularized equations (0.1) says that the ki¬ 
netic energy variation in a given volume element A of the fluid, in a time 
interval [to,ti\, plus the energy dissipated therein by friction, equals the 
sum of the kinetic energy that in the time interval t £ [to, t{\ enters in the 
volume element plus the work performed by the pressure forces (on the 
boundary element) plus the work of the volume forces (none in our case). 
The analytic form of this relation is simply obtained by multiplying both 
sides of the first of the (0.1) by u and integrating on the volume element 
A and over the time interval [to,ti\. 

The relation that one gets can be generalized to the case in which the vol¬ 
ume element has a time dependent shape. And an even more general re¬ 
lation can be obtained by multiplying both sides of (0.1) by (f(x, t)u(x , t) 
where if is a x (0, s]) function with f(x,t) zero for t near 0 (here 

s is a positive parameter). 

Energy conservation in a sharply defined volume A and time interval 
t G [to,C] can be obtained as limiting case of choices of f in the limit 
in which it becomes the characteristic function of the space-time volume 
element A x [to,ti\. 
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1.3 


1.4 


Making use of a regular function (f(x, t ) is useful, particularly in the rather 
“desperate” situation in which we are when using the theory of Leray. 
The “solutions” u (obtained by removing, in (0.1), the regularization, 
i.e. letting A —> oo) are only “weak solutions”. Therefore, the relations 
that are obtained in the limit Acjoc can be interpreted as valid only 
after suitable integrations by parts that allow us to avoid introducing 
derivatives of u (whose existence is not guaranteed by the theory) at the 
“expense” of differentiating the “test function” (f. 

Performing analytically the computation of the energy balance, described 
above in words, in the case of the regularized equation (0.1) and via a 
few integrations by parts 2 we get the following relation 


\ [ d£\u(^s)\ 2 (p(£,s) + v [ dt f (f(^t)\du(^t)\ 2 dx = 
z Jn Jo Jn 


■((ft + v Acp)\u\ 2 + -\u\ 2 (u)x ■ dtp + pu ■ dcp dtdl- 



o Jn 


( 1 . 1 ) 


where (f t =dt<f and u = u x is in fact depending also on the regularization 
parameter A; here p is the pressure p = — A ~ 1 didj(uiUj). 

Suppose that the solution of (0.1) with fixed initial datum u 0 converges 
(weakly in L 2 ), for A —» 00 , to a “Leray solution” u possibly only over a 
subsequence A n —> oc. 

The (1.1) implies that any (in case of non uniqueness) Leray solution u 
verifies the energy inequality: 


'n 


\u(i,s)\ 2 (f(£,s)d£ + is / (f(^t)\du(^t)\ 2 d^dt < 

Jt<s j 


< 


h < 8 Jn 


1 


1 


( 1 . 2 ) 


-((ft + v A(f)\u\ z + -1 u\u ■ d(f + pu ■ d(f d£ dt 


where the pressure p is given by p = — A 1 didj(uiUj) = 

—A ~ 1 dd(uu) = —A ~ 1 (du) 2 - 


2 


The solutions of (0.1) are C°°(Q X [0, 00 )) so that there is no need to justify integrating 
by parts. 
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Remarks: (1) It is important to remark that in this relation one might 
expect the equality sign: as we shall see the fact that we cannot do better 
than just obtaining an inequality means that the limit necessary to reach 
a Leray solution can introduce a “spurious dissipation” that we are simply 
unable to understand on the basis of what we know (today) about the 
Leray solutions. 

(2) The above “strange” phenomenon reflects our inability to develop a 
complete theory of the Navier-Stokes equation, but one can conjecture 
that no other dissipation can take place and that a (yet to come) com¬ 
plete theory of the equations could show this. Hence we should take the 
inequality sign in (1.2) as one more manifestation of the inadequacy of 
the Leray’s solution. 

The proof of (1.2) and of the other inequalities that we shall quote and 
use in this section is elementary and based, c./.r. problem [15] below, on 
a few general “kinematic inequalities” that we now list (all of them will 
be used in the following but only (S) and (CZ) are needed to check (1.2)). 

§2. Kinematic inequalities. 

A first “kinematic” inequality, i.e. the first inequality that we shall need 
and that holds for any function /, is 3 

(P) Poincare inequality: 


f dx\f - F\ a < C? r 3 ~ 2 “ ( f 

B r \ ^3j- 


a 


dx \df\ 


1 < « < 2 f 2 - 1 ) 


where F is the average of / on the ball B r with radius r and is a 
suitable constant. We shall denote (2.1) by (P). 


A second kinematic inequality that we shall use is 


The inequalities should be regarded as inequalities for C°° functions; they can be 
extended to the appropriate Sobolev spaces by continuity. 
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(S) Sobolev inequality: 


2.2 



u\ q dx < Cq 





q/2—a 

+ 


_|_ r 


—2a 



y/ 2 , 

\u\ 2 dx ) 


J 


if 2 < q < 6, 


“= 


( 2 . 2 ) 


where B r is a ball of radius r and the integrals are performed with respect 
to dx. The C^ is a suitable constant; the second term of the right hand 
side can be omitted if u has zero average over B r . We shall denote (2.2) 
by (S), [So63], 

A third necessary kinematic inequality will be 
(CZ) Calderon-Zygmund inequality: 


2.3 




1 didj)(uiUj)\ q d£ t < Cq / \u\ 2q dt; , 1 < q < oc (2.3) 


which we shall denote (CZ): here Si is the torus of side L and C^ is a 
suitable constant, [St93]. 

And finally 

(H) Holder inequality: 


2.4 


/l/2 • • • fr 


2—1 


1/. 


1 n 


E 

2—1 


= 1 


Vi 


(2.4) 


which we shall denote (H): the integrals are performed over an arbitrary 
domain with respect to an arbitrary measure (of course the same for all 
integrals). 

Remark: The (H) are a (simple) extension of the Schwartz-Holder in¬ 
equalities; the (P) is a simple inequality for a = 1, c./.r. problem [13]; 
while (P), (S) (mainly in the cases a = | or q = 6) and (CZ) are less 
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elementary and we refer to the literature, footnote at p.43 [So63], [St93], 
[LL01] for their proofs. 

An important consequence of the inequalities is 

Proposition 1: Let u be a Leray solution verifying (therefore) the a 

T 

priori bounds in (0.2): J ^ \u(x, t)\ 2 dx < Eq and f Q dt f ^ \du(x, t)\ 2 dx < 
Eqis~ 1 then 

f dt f dx\u\ 10 / 3 + f dt f dx\p\ 5 / 3 < Co~ 1 Eq ^ 3 (2.5) 
Jo Jn Jo Jn 

where C can be chosen Cfo (1 + C%). 

3 3 

proof: Apply (S) with q = ^ and a = 1: 


u\ 3 dx < Cfo ( / (du) 2 dx 


u 2 dx 


<CfoE* / \du\ 


( 2 . 6 ) 

hence integrating over t between 0 and T using also the second a priori 
estimate, we find 


dt / \u\~dx < Cfo Eq / dt dx(du) 2 <Cfor , ~ 1 EQ +3 


(2.7) 


while the (CZ) yields: dx \p\ 3 < Cf dx \u\ 3 ° which, integrated over 

t and combined with (2.7), gives the announced result. 

§3. Pseudo Navier Stokes velocity—pressure pairs. Scaling op¬ 
erators. 


As already mentioned the CKN theory will not fully use that u verifies 
the Navier-Stokes equation: in order to better realize this (unpleasant) 
property it is convenient to define separately the only properties of the 
Leray solutions that are really needed to develop the theory, i. e. to obtain 
an estimate of the fractal dimension of the space-time singularities set 
Sq. This leads to the following notion 
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Definition a: (pseudo NS velocity field): Lett —» (u(-, t), p(-, t)) be a 
function with values in the space of zero average square integrable “veloc¬ 
ity” and “pressure” fields on fh Suppose that for each cp G C°°(Llx (0,71) 
with ip(x,t) vanishing for t near zero the following properties hold. For 
each T < oc and s < T: 

(a) / udx = 0, d ■ u = 0, P = — / djdjA -1 (ujUj) 

Jn id 

( b ) f dt f dx\u \ 10 / 3 + f dt f dx\p \ 5 ^ 3 < oo 

Jo Jn Jo Jn 

(c) i I dx 

<-j 

J t<s 

(3.1) 

Then we shall say that the pair (u,p) is a pseudo NS velocity and pressure 
pair. The singularity set in the time interval [0, T] of (u,p) will be defined 
as the set S o of the points (x,t) G Tl x [0, T] that do not admit a vicinity 
U where \u\ is bounded . 4 

The name given to the set So is justified by a general result on the theory 
of NS equations which shows that if a Leray’s solution of the NS equations 
is essentially bounded in a neighborhood of a space time point then it is 
C°° near such point. 

Proposition 2: (velocity is unbounded near singularities) : Letu(x,t) be 
a Leray’s solution of the NS equation in L- 2 . Given to > 0 suppose that 
\u(x,t)\ < M, (x,t) G U p { xq , to ) = sphere of radius p (p < to) around 
(xqUo), for some M < oo: then u G C°°(U p / 2 (x 0 ,to)). 

Remarks: 


\u(x, s)\ 2 (p(x, s) + v / / ip{x, t)\du\ 2 dxdt < 

Jt<s JQ 


'n 


-((ft + o Acp)\u \ 2 H— \o\ 2 o • d<T Tpu -dp 


dx dt 


Here we mean bounded outside a set of zero measure in U or, as one says, essentially 
bounded because it is clear that, being u,p in they are defined up to a set of 

zero measure and it would not make sense to ask that they are bounded everywhere 
without specifying which realization of the functions we take. 
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(1) This means that the only way a singularity can manifest itself, in a 
Leray solution of the NS equations, is through a divergence of the velocity 
held itself. For instance it is impossible to have a singular derivative hav¬ 
ing the velocity itself unbounded. Hence, if d > 3 velocity discontinuities 
are impossible (and even less so shock waves), for instance. Naturally if 
u(x, t ) is modified on a set of points (x, t) with zero measure it remains 
a weak solution (because the Fourier transform, in terms of which the 
notion of weak solution is defined, does not change), hence the condition 
| u(x, t) | < M for each (x, t) G [/y(x 0 , t 0 ) can be replaced by the condition: 
for almost all (x, t) G U\(x 0 ,to). 

(2) The above result is not strong enough to overcome the difficulties 
of a local theory of regularity of the Leray weak solutions. Therefore 
one looks for other results of the same type and it would be desirable 
to have results concerning regularity implied by a priori informations on 
the vorticity. We have already seen that bounded total vorticity implies 
regularity: however it is very difficult to go really beyond; hence it is 
interesting to note that also other properties of the vorticity may imply 
regularity. A striking result in this direction, although insufficient for 
concluding regularity (if true at all) of Leray weak solutions, is in [CF93]. 

(3) For a proof of the above (Serrin’s) theorem see [Ga02], proposition 
IV in section 3.3. 

The remaining part of this section will concern the general properties 
of the pseudo NS pairs and their regularity at a given point (x, t): it 
will not have more to do with the velocity and pressure fields that solve 
the Navier-Stokes equations. It is indeed easy to convince oneself that 
the (3.1), in spite of the arbitrariness of (p, are not equivalent , not even 
formally, to the Navier-Stokes equations, and they pose far less severe 
on u,p restrictions. We should not be surprised, therefore, if it turned 
out possible to exhibit pseudo NS pairs that really have singularities on 
“large sets” of space-time. In a way it is already surprising that the 
pseudo NS fields verify the regularity properties discussed below. 
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3.2 


3.3 


3.4 


3.5 


The analysis of the latter properties (of pseudo NS fields) is based on the 
mutual relations between certain quantities that we shall call “ dimen¬ 
sionless operators ” relative to the space-time point Qr 0 ,£o) 

Definition b: (dimensionless “operators” for NS) Let (x 0 , to ) G fi x 
(0, oo) and consider the sets 5 

A r (t 0 ) ={*| 1 1 - t 0 1 < r 2 ^ 1 } 

B r {x o) ={£l l£-£ol <r} = B r 

Qr(x 0 ,t 0 ) ={{£,$) | g-^ol < r, \tf - to | < r 2 v~ 1 } 

Qr{x_ Q, to) — N r (t o) x B r (x_ q) = Q r 

define: 

(i) “dimensionless kinetic energy operator” on scale r: 

1 


(3.2) 


A ( r ) = — sup / HU)\ df 

^ ^ \t — to|<n — 1 r 2 JB r 


(3.3) 


and we say that the dimension of A is 1 : this refers to the factor r ~ 1 
that is used to make the integral dimensionless. 

(ii) “local Reynolds number” averaged on scale r: 


5(r) = — f dddf\du 


(3.4) 


and we say that the dimension of 6 is 1 : this refers in general to the 
power —a to which r has to be raised so that an expression becomes di¬ 
mensionless: in this case a = 1. 

(Hi) “dimensionless energy flux” on scale r: 


G(r) = 


1 


v 2 r 2 J Qt 


dddf u 


(3.5) 


5 If r > L/2 this is interpreted as B r — O. If r 2 v 1 > to then A r (to) is interpreted as 
0 < t < to + r 2 u~ i . 
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The dimension of G is 2. 

(iv) “dimensionless pressure power” forces on scale r: 

J(r) = 2 2 f dfdd |a| \p\ (3-6) 

«/ Q r 

The dimension of J is 2. 

(v) “dimensionless non locality” on scale r: 

r ~ 13 / 4 f / f \ 5/4 

K ( r ) = 3/2 / ( 3 ' 7 ) 

The dimension of K is 13/4. 

(vi) “dimensionless intensity ” on scale r: 


3.8 


3.9 


S'(r) = i,- 7 /3 r -5/3 f (| m | 10/3 + | p |5/3^^ (3.8) 

J Q r 


where the pressure is always defined by the expression p = 
> 2 ,; (uiUj). The dimension of S is 5/3. 

Remarks: (1) The 2 L(r),... are not operators in the common sense of 
functional analysis. Their name is due to their analogy with the quan¬ 
tities that appear in problems that are studied with the methods of the 
“renormalization group” (which, also, are not operators in the common 
sense of the words). Perhaps a more appropriate name could be “di¬ 
mensionless observables”: but we shall call them operators to stress the 
analogy of what follows with the methods of the renormalization group. 

( 2 ) The A{r), G{r ), J(r), K(r ), S{r) are in fact estimates of the quantities 
that their name evokes. We omit the qualifier “estimate” when referring 
to them for brevity. 

(3) The interest of (i)q-(iv) becomes manifest if we note that the en¬ 
ergy inequality (3.1) can be expressed in terms of such quantities if cp is 
suitably chosen. Indeed let 


V = X{x,t) 


exp — ( 


_ (T g 0 ) _\ 

4(i/po— t)+2r 2 ) ) 


(A.Tru{t — to) + 87 rr 2 ) 3 / 2 


(3.9) 
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3.10 


3.11 


where yQr, t) is C 00 and has value 1 if (x,t) E Q r /2 and 0 if (x,t) Q r . 
Then there exists a constant C > 0 such that 


11 C 10 1 C 

M < —, \d<p\ < —, 

ry» u ry‘ j 

1^1 > 


\dt(f + \ < _ 1 r , everywhere 

z/ 1 r t> 

if (x, t ) E Q r / 2 

(3.10) 


Hence (3.1) implies 



(3.11) 


and, since Jg \u \ 2 < C (Jq \u\ 3 ) 2 ^ 3 (o V 5 ) 1 / 3 with a suitable C, it fol¬ 
lows that for some C 


n(§) + <5(5) < C (G(r) 2 / 3 + G(r) + J(r)) (3.12) 

(4) Note that the operator 5(r) is an average of the “local Reynolds’ 
number” r f A df\du\ 2 . 

(5) The operator (v) appears if one tries to bound J(|) in terms of 
A{r) + 5(r): such an estimate is indeed possible and it will lead to the 
local Scheffer theorem discussed in the next section. 

§4. The theorems of Scheffer and of Caffarelli—Kohn—Nirenberg. 

We can state the strongest results known (in general and to date) about 
the regularity of the weak solutions of Navier Stokes equations (which 
however hold also for the pseudo Navier Stokes velocity-pressure pairs). 

Theorem I: (upper bound on the dimension of the sporadic set of singular 
times for NS, ( “Scheffer’s theorem”)): There are two constants e s ,C > 0 
such that if G(r) + J(r) + K (r) < e s for a certain value of r, then u is 
bounded in Qq(x 0 ,to) : 

\u(x,t)\ < C -, (x,t) E Q^{x 0 ,to), almost everywhere (4.1) 
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having set v = 1. 

Remarks: (1) c.f.r. problems [5] 4-[11] for a guide to the proof. 

(2) This theorem can be conveniently combined, for the purpose of 
checking its hypotheses, with the inequality: J(r) + G{r) + K (r) < 
C (yS (r) 9 / 10 + S'(r) 3 / 4 ), which follows immediately from inequality (H) 
and from the definitions of the operators, with a suitable C. 

(3) In other words if the operator S(r ) is small enough then (t 0 No) is a 
regular point. 

(4) This will imply that the fractal dimension of the space-time singu¬ 
larities set is < 5/3. In fact, see section 5 below, an a priori estimate on 
the global value of an operator with dimension a implies that the Haus¬ 
dorff ’ measure of the set of points around which the operator is large 
does not exceed a; here the operator S(r ) has dimension 5/3 and there¬ 
fore together with the a priori bound (2.5) it yields and estimate <5/3 
for the Hausdorff dimension of the singularity set. This also justifies the 
introduction of the operator S(r). 

It is easy, in terms of the just defined operators, to illustrate the strat¬ 
egy of the proof of the following theorem which will immediately imply, 
via a classical argument reproduced in section 5 below, that the fractal 
dimension of the space time singularities set So for a pseudo NS held is 
< 1 and that its 1-measure of Hausdorff /ii(S'o) vanishes. 

Theorem II: (sufficient condition for local regularity space-time (“CKN 
theorem”)) There is a constant e c r, n such that if (u,p) is a pseudo NS 
pair of velocity and pressure fields and 

, t')\ 2 dx'dt' = limsup 5(r) < £ c k n (4-2) 

r—>0 

then u(x',t'), p(x',t') are C°° in the vicinity of ( XqRo )- 6 

6 This means that near ( x , t ) the functions u(x', t'),p(x ', t') coincide with C°° functions 
apart form a set of zero measure (recall that the pseudo NS fields are defined as fields 
in L 2 CI)). 


lim sup — / | du(x' 

r^o vr jQ r (x 0 ,t 0 ) ' 
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For fixed (x 0 ,to), consider the “sequence of length scales”: r n = L2 n , 

with n = 0,-1, —2,_ We shall set a n = A{r n ), K n = K^ 5 , j n = J n , 

g n = G^/ 3 , = 4(r n ) which is a natural definition as it will shortly 

appear. And define X n = (ct n , n n , j n , g n ) £ R\. Then the proof of 
this theorem is based on a bound that allows us to estimate the size of 
X_ n . dehned as the sum of its components, in terms of the size of X_ n+p 
provided the Reynolds number S n+P on scale n + p is < 5. 

The inequality will have the form (if p > 0 and 0 < 5 < 1) 

X n <B p (X n+p ;5) (4.3) 

where B p f; 5) is a map of the whole into itself and the inequality has to 
be understood “component wise”, i. e. in the sense that each component 
of the l.h.s. is bounded by the corresponding component of the r.h.s. We 
call |X| the sum of the components of X £ R\. 

The map B p {-] 6 ), which to some readers will appear as strongly related to 
the “ beta function ” for the “running couplings” of the “renormalization 
group approaches”, 7 will enjoy the following property 

Proposition 3: Suppose that p is large enough; given p > 0 there exists 
6 p (p) > 0 such that if 5 < S p (p ) then the iterates of the map B p (-;5) 
contract any given ball in R\, within a finite number of iterations, into 
the ball of radius p: i.e. \B p (X_]5)\ < p for all large k’s. 

Assuming the above proposition the main theorem II follows: 

proof: Let p = e s , c.f.r. theorem I, and let p be so large that the above 
proposition holds. We set £ c k n = 5 P {£ S ) and it will be, by the assumption 
(4.2), that S n < £ c k n fo r all n < no for a suitable no (recall that the scale 
labels n are negative). 


Indeed it relates properties of operators on a scale to those on a different scale. Note, 
however, that the couplings on scale n, i.e. the components of X , provide information 
on those of X n +p rather than on those of X_ n p as usual in the renormalization group 

methods, see [BG95]. 
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Therefore it follows that \B p (X nQ \ £ r .k.n)\ < £ s for some k. Therefore by 
the theorem I we conclude that (x 0 ,to) is a regularity point. 

Proof that the renormalization map contracts. 

Proposition 3 follows immediately from the following general “Sobolev 
inequalities” 

(1) “Kinematic inequalities”: i.e. inequalities depending only on the 
fact that u is a divergence zero, average zero and is in Z/2(f4) and 
p = —A ~ 1 (du) 2 

J n <C (2 - p/5 A 1 £ p G 1 7 / 5 K*!? p + 2 2 ”A 1 J 2 p S n+p ) 

K n < C (2 ~” /2 K n+p + (4.4) 

G 2 J Z < G (2- 2 M„ +p + 2 2 My + 2 p O 

where C denotes a suitable constant ( independent on the particular 
pseudo NS field). The proof of the inequalities (4.4) is not difficult, as¬ 
suming the (S,H,CZ,P) inequalities above, and it is illustrated in the 
problems [1], [2], [3]. 

(2) “ Dynamical inequality”: i.e. an inequality based on the energy in¬ 
equality (c) in (3.1) which implies, quite easily, the following u dynamic 
inequality ” 8 

A n < C (2 p Gj +p + 2 p A n + p 5 n + p + 2 p J n + p ) (4-5) 

whose proof is illustrated in problem [4], 

proof of proposition: Assume the above inequalities (4.4), (4.5) and set¬ 
ting a n = A n ,K n = I<n /5 Jn = Jn; Qn = G % 3 , S n+P = 6 and, as above, 
A’ — {cx n , tz n , jn , Qn)• The r.h.s. of the inequalities defines the map 
B P (X-S). 


We call it “dynamic” because it follows from the energy inequality, i.e. from the 
equations of motion. 
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If one stares long enough at them one realizes that the contraction prop¬ 
erty of the proposition is an immediate consequence of 

(1) The exponents to which £ = 2~ p is raised in the various terms are 
either positive or not; in the latter cases the inverse power of £ is always 
appearing multiplied by a power of 5 n + p which we can take so small to 
compensate for the size of £ to any negative power, except in the one 
case corresponding to the last term in (4.5) where we see e~ 1 without any 
compensating 5 n + p . 

(2) Furthermore the sum of the powers of the components of X_ n in each 
term of the inequalities is always < 1: this means that the inequalities 
are “almost linear” and a linear map that “bounds” B p exists and it is 
described by a matrix with small entries except one off-diagonal element. 
The iterates of the matrix therefore contract unless the large matrix ele¬ 
ment “ill placed” in the matrix: and one easily sees that it is not. 

A formal argument can be devised in many ways: we present one in which 
several choices appear that are quite arbitrary and that the reader can 
replace with alternatives. In a way one should really try to see why a 
formal argument is not necessary. 

The relation (4.5) can be “iterated” by using the expressions (4.4) for 
J n + P and then the first of (4.4) to express G\[+ p in terms of A n+ 2 P 
with n replaced by n + p: 


a n < C (2 p Oi n+ 2 P + 2 3p dJ + 2 P a r!+2 P + 

+ 2 p / 5 (a n+ 2pKn+2p) 1 / 2 + 2 7p / 5 4 n+ 2p<a^+2p K: i+2p+ 

T 2 ^d n j r 2p(y,ri-\-2p) 

It is convenient to take advantage of the simple inequalities ( ab )^ < 
za + z~ 1 b and a x < 1 + a for a, 6, z, x > 0, x < 1. 

The (4.6) can be turned into a relation between a n and a n + p , n n + P by re¬ 
placing p by \p . Furthermore, in the relation between a n and a n+pi K n+P 
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obtained after the latter replacement, we choose z = 2 p / 5 to disentangle 
2 p / 10 (a n+p K, n+p ) 1 / 2 we obtain recurrent (generous) estimates for a n , k 2 


4.7 


«n < C (2 - p / 10 a„ +p + 2 3p / 1 V +P + £° +p ) 

«n < C (2 -4p / 5 K n+p + £ +p ) 


£n+p — 2 3p (5 n -(- p (a n + p + K n + P + 1) 

C+p = 2 3p <5„ +p a„ +p 


(4.7) 


We fix p once and for all such that 2 p / 10 C < 77 . 

Then if C2 3p S n is small enough, i.e. if 5 n is small enough, 

say for 5 n < 5 for all \n\ > n , the matrix M = 

f 2~ p / 10 + 2 3p 5 n+p 2 3p / 10 + 2 3p 6 n+p \ ... . + . . 

Cl 0 2 ^ 4 p / 5 + 2 3p S J W1 ^ ^ iave ^ w0 ei § enva f ues 

< and iteration of (4.6) will contract any ball in the plane a, k to the 

ball of radius 2 5. 

If a n , K n are bounded by a constant 5 for all \n\ large enough the (4.4) 
show that also g n ,jn are going to be eventually bounded proportionally 
to 5. 

Hence by imposing that 5 is so small that \X n \ = a n + K n + jn + 9n < p 
p we see that proposition 3 holds (hence theorem 2 as a consequence of 
theorem 1 ). 


§5. Fractal dimension of singularities of the Navier—Stokes equa¬ 
tion, d = 3. 


Here we ask which could be the structure of the possible set of the sin¬ 
gularity points of the solutions of the Navier-Stokes equation in d = 3. 
The answer is an immediate consequence of theorem II and we describe 
it here for completeness: the technique is a classic method (Almgren) to 
link a priori estimates to fractal dimension estimates. 

It has been shown already by Leray that the set of times at which a 
singularity is possible has zero measure (on the time axis), see §3.4 in 
[Ga02], 
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Obviously sets of zero measure can be quite structured and even large in 
other senses. One can think to the example of the Cantor set which is non 
denumerable and obtained from an interval I by deleting an open con¬ 
centric subinterval of length 1/3 that of I and then repeating recursively 
this operation on each of the remaining intervals (called n-th generation 
intervals after n steps); or one can think to the set of rational points 
which is dense. 

(A) Dimension and measure of Hausdorff. 

An interesting geometric characteristic of the size of a set is given by the 
Hausdorff dimension and by the Hausdorff measure, c.f.r. [DS60], p.174. 

Definition c: (Hausdorff a-measure): The Hausdorff a-measure of a 
set A contained in a metric space M is defined by considering for each 
5 > 0 all coverings C§ of A by closed sets F with diameter 0 < d{F ) < d 
and setting 

Tcx(A) = lim inf ^ d(F) a (5.1) 

5 FeC s 

Remarks: (1) The limit over 6 exists because the quantity inf^ ... is 
monotonic nondecreasing. 

(2) It is possible to show that the function defined on the sets A of 
M by A — > n a {A) is completely additive on the smallest family of sets 
containing all closed sets and invariant with respect to the operations of 
complementation and countable union (which is called the a-algebra E 
of the Borel sets of M), c.f.r. [DS60]. 

One checks immediately that given A £ E there is a c such that 

H a (A ) = oo if a < a c , p a (A) = 0 if a > a c (5-2) 

and it is therefore natural to set up the following definition 

Definition d: (Hausdorff measure and Hausdorff dimension): Given a 
set A C R d the quantity a c , (5.2), is called Hausdorff dimension of A, 
while p ac (A) is called Hausdorff measure of A. 
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It is not difficult to check that 

(1) Denumerable sets in [0,1] have zero Hausdorff dimension and measure. 

(2) Hausdorff dimension of n-dimensional regular surfaces in R d is n and, 
furthermore, the Hausdorff measure of their Borel subsets defines on the 
surface a measure p ac that is equivalent to the area measure p: namely 
there is a p{x ) such that p ac (dx) = p{x ) p(dx). 

(3) The Cantor set, defined also as the set of all numbers in [0,1] which 
in the representation in base 3 do not contain the digit 1, has 

a c = log 3 2 (5.3) 


as Hausdorff dimension. 9 

(B) Hausdorff dimension of singular times in the Navier-Stokes solutions 
(d = 3). 

We now attempt to estimate the Hausdorff dimension of the sets of 
times t < T < oo at which appear singularities of a given weak solu¬ 
tion of Leray, i.e. a solution of the type discussed in (0.1). Here T is an 
a priori arbitrarily prefixed time. 

We need a key property of Leray’s solutions, namely that if at time to 
it is Ji(t 0 ) = L~ l f(du) 2 dx < oo, i.e. if the Reynolds number R(to) = 

d j~ 

Mto) 1/2 /V c = V 1 /V c with V c = vL- 1 is < Too, then the solution stays 


Indeed with 2 n disjoint segments with size 3 n , uniquely determined (the n —th gen¬ 
eration segments), one covers the whole set C: hence 

Ua,S d — inf y d(F) a <1 if a — ao = log 3 2 
c s 

Fec s 

and g ao {C) <1: i.e. /r a (C) = 0 if a > ao- Furthermore, c.f.r. problem [16] below, if 
a < ao one checks that the covering C° realizing the smallest value of d(F) a 

with 5 — 3~ n is precisely the just considered one consisting in the 2 n intervals of 
length 3 -n of the n—th generation and the value of the sum on such covering diverges 
for n —> oo. Hence /i Q -(C') = oo if a < ao so that ao = a c and p ac (C) = I- 
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5.4 


regular in a time interval (to, to + r] with (see proposition II in §3.3 of 
[Ga02], eq. (3.3.34)): 


T 

t = F --— 

R(toY' 



V 


(5.4) 


From this it will follow, see below, that there are A > 0,7 > 0 such that 
if 

liminf f^Y f —R 2 (id) < A (5.5) 

cr— \T C J J t _ a (j 

then t > a and the solution is regular in an interval that contains t so 
that the instant t is an instant at which the solution is regular. Here, as 
in the following, we could fix 7 = 1 / 2 : but 7 is left arbitrary in order to 
make clearer why the choice 7 = 1/2 is the “best”. 

We first show that, indeed, from (5.5) we deduce the existence of a se¬ 
quence (i t 0 such that 


5.6 



did 






(5.6) 


therefore, the l.h.s. being a time average, there must exist do i E (t — <Ji,t) 
such that 

^ 0 ,)<.4(fr < 5 - 7 ) 

and then the solution is regular in the interval (do i, do i + rf) with length 
Ti at least 

t? T (Vi/Tc) 21 ^ 

Ti = FT C - — - > ai (5.8) 

provided 7 < 1 / 2 , and is small enough and if A is small enough (if 
7 =| then this means 2A 2 < F). Under these conditions the size of the 
regularity interval is longer than cq and therefore it contains t itself. 

It follows that, if t is in the set S of the times at which a singularity is 
present, it must be 

liminf [ —R 2 (d)>A if t £ S (5.9) 

(t^o \T C J J t _ G a 
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5.10 


5.11 


5.12 


i.e. every singularity point is covered by a family of infinitely many in¬ 
tervals F with diameters a arbitrarily small and satisfying 



didR 2 {d) > 




(5.10) 


From Vitali’s covering theorem (c./.r. problem [19]) it follows that, given 
5 > 0 , one can find a denumerable family of intervals F 2 ,..with 
Fi = (ti — (Ti,ti ), pairwise disjoint and verifying the (5.10) and 0 { < 5/ 4, 

def 

such that the intervals 5 Fi = (^ — 7<t*/ 2, A + 5cp,/2) (obtained by dilating 
the intervals iy by a factor 5 about their center) cover S 


S C Uf 5 Fi 


(5.11) 


Consider therefore the covering C generated by the sets 5 Fi and compute 
the sum in (5.1) with a = 1 — 7 : 



2 5 1 - 7 ^ 

< avt c G 


r 0151-7 


*T 


did R 2 (d ) < 00 


(5.12) 


where we have made use of the a priori estimates on vorticity ( 0 . 2 ) and 
we must recall that 7 < 1/2 is a necessary condition in order that what 
has been derived be valid (c./.r. comment to (5.8)). 

Hence it is clear that for each a > 1/2 it is p a {S) < oc (pick, in fact, 
a = 1 — 7 , with 7 < 1 / 2 ) hence the Hausdorff dimension of S is a c < 
1 / 2 . Obviously the choice that gives the best regularity result (with the 
informations that we gathered) is precisely 7 = 1 / 2 . 

Moreover one can check that H\/ 2 {S) = 0: indeed we know that S has 
zero measure, hence there is an open set G D S with measure smaller 
than a prefixed e. And we can choose the intervals F t considered above 
so that they also verify F t c G: hence we can replace the integral in 
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the right hand side of (5.12) with the integral over G hence, since the 
integrand is summable, we shall find that the value of the integral can be 
supposed as small as wished, so that fii/ 2 (S) = 0. 

(C) Hausdorff dimension in space-time of the solutions of NS, (d = 3). 

The problem of which is the Hausdorff dimension of the points (x , t ) G 
H x [0, T] which are singularity points for the Leray’s solutions is quite 
different. 

Indeed, a priori , it could even happen that, at one of the times t G S 
where the solution has a singularity as a function of time, all points 
( x,t ), with igO, are singularity points and therefore the set So of the 
singularity points thought of as a set in space-time could have dimension 
3 (and perhaps even 3.5 if we take into account the dimension of the 
singular times discussed in (B) above). 

From theorem II we know that if 


rt+r 2 /2u n in 

limsupr -1 / / — df{d_u) 2 <e (5.13) 

0 J t—r 2 /2v JS(x_,r ) ^ 

then regularity at the point ( x , t ) follows. 

It follows that the set So of the singularity points in space-time can be 
covered by sets C r = S(x,r) x (t — r 2 u~ l ,t + with r arbitrarily 

small and such that 


5.14 


1 

rv 



dd 


S(x,r) 


dx ( du) 2 > £ 


(5.14) 


which is the negation of the property in (5.13). 

Again by a covering theorem of Vitali (c./.r. problems [16],[17]), we can 

2 2 

find a family Fi of sets of the form Fi = S(x i: ri ) x (A — 
pairwise disjoint and such that the sets 6iq= set of points (x',t') at 
distance < 6r* from the points of F, h covers the singularity set Aq. 10 One 


10 


Here the constant 5, as well as the other numerical constants that we meet below like 
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5.15 


can then estimate the sum in (5.1) for such a covering, by using that the 
sets F % are pairwise disjoint and that 5 Fi has diameter, if maxrj is small 
enough, not larger than 18rp 


y~Vl8rj) < —f (du) 2 dfdt < — f f (du) 2 dfdt < oo (5.15) 


■>T 


i.e. the 1-measure of Hausdorff /ii(So) would be < oo hence the Hausdorff 
dimension of So would be < 1. 

Since So has zero measure, being contained in Si x S where S is the set of 
times at which a singularity occurs somewhere, see (5.9), it follows (still 
from the covering theorems) that in fact it is possible to choose the sets 
Fi so that their union U is contained into an open set G which differs 
from So by a set of measure that exceeds by as little as desired that of So 
(which is zero); one follows the same method used above in the analysis 
of the time-singularity. Hence we can replace the last integral in (5.15) 
with an integral extended to the union U of the F's: the latter integral 
can be made as small as wished by letting the measure of G to 0. It 
follows that not only the Hausdorff dimension of So is < 1, but also the 
Mi (S 0 ) = 0. 


Remarks: (1) In this way we exclude that the set Sq of the space-time 
singularities contains a regular curve: singularities, if existent , cannot 
move along trajectories (like flow lines) otherwise the dimension of S 
would be 1 > 1/2) nor they can distribute, at fixed time, along lines and, 
hence, in a sense they must appear isolates and immediately disappear 
(always assuming their real existence). 

(2) A conjecture (much debated and that I favor) that is behind all our 
discussions is that if the initial datum u° is in C°°(Sl) then there exists 
a solution to the Navier Stokes equation that is of class C°° in (x, t)”, 
i.e. So = 0! The problem is, still, open: counterexamples to the conjec¬ 
ture are not known (i.e. singular Leray’s solutions with initial data and 


5, 6, 18 have no importance for our purposes and are just simple constants for which 
the estimates work. 
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external force of class C°°) but the matter is much debated and different 
alternative conjectures are possible (c./.r. [PS87]). 

(3) In this respect one should keep in mind that if d > 4 it is possible 
to show that not all smooth initial data evolve into regular solutions: 
counterexamples to smoothness can indeed be constructed, c./.r. [Sc77]). 


Problems. The dimensional bounds of the CKN theory. 

In the following problems we shall set v — 1, with no loss of generality, thus fixing 
the units so that time is a square length. The symbols (u, p) will denote a pseudo 
NS field, according to definition 1. Moreover, for notational simplicity, we shall set 
A p = A(p), G p = G(p),..., and sometimes we shall write A rn ,Gr n ... as A n ,... 
with an abuse that should not generate ambiguities. The validity of the (3.1) for 
Leray’s solution is checked in problem [15], at the end of the problems section, to 
stress that the theorems of Scheffer and CKN concern pseudo NS velocity—pressure 
fields: however it is independent of the first 14 problems. There will many constants 
that we generically denote C: they are not the same but one should think that they 
increase monotnically by a finite amount at each inequality. The integration elements 
like dx and dt are often omitted to simplify the notations and they should be easily 
understood from the integration domains. 

[1]: Let p = r n _|_p and r = r n , with r n — L2 n , c./.r. lines following (4.2), and apply 
(S),(2.2), with q = 3 and a — to the field u, at t fixed in A r and using definition 2 
deduce 



Infer from the above the third of (4.4). (Idea: Let \u\ 2 be the average of y? on the 
ball Bp ; apply the inequality (P), with a = 1, to show that there is C > 0 such that 


dx I ur < 


'B, 


dx 


Br 


m - Mp 


+ \u 


—I p 


dx < 


B r 


B 


<Cp / dx |u| \du | + C ( — j / dx \u\ 2 < Cp^^A 1 ^ 2 ( / dx\du\ 2 ) + 


Br 


B r 


1/2 


+c w pAp 

where the dependence from t E A r is not explicitly indicated; hence 

j dx\u\ 3 <C(rp~ 1 ) 3 A^ 2 + C(p |3//4 T p 9//4 r _ J,,/2 )Ap^ 4 ^ f dx\du\ 2 ^j 


3/4 
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then integrate both sides with respect to t £ A r and apply (H) and definition 2.) 

[ 2 ]: Let p < 1 be a non negative C°° function with value 1 if \x\ < 3p/4 and 0 if 
|s| > 4p/5; we suppose that it has the “scaling” form p = <p\(x/p with pi > 0 a C°° 
function fixed once and for all. Let B p be the ball centered at x with radius p; and 
note that, if p = r n _|_ p and r = r n , pressure can be written, at each time (without 
explicitly exhibiting the time dependence), as p(x) — p'(x) + p"(x) with 


p'(x) = 


47r 


B r 


7—-—tp(p) A(p(p) dy + 

|x —y| — — 27 t 




x-y 
\x — y | £ 


9<p(y)p(y) dy 


p" (x ) =— 
v - ; 4 tt 


l B \x-y\ 


■p(y) (du(y)) ■ ( du(y))dy 


if |x| < 3p/4; and also |p 7 (x)| < Cp 3 dy\p(y)\ and all functions are evaluated 

J ±3 p — — 

at a fixed t £ A r . Deduce from this remark the first of the (4.4). (Idea: First note 
the identity p = — (47 t)~ ] j b \x_ — y\~ ] A (pp) for x £ B r because if x £ 5 3p / 4 it 

is <pp = p. Then note the identity A (<pp) = p A <p + 29p • 9<p + <p Ap and since 
Ap — —d ■ (u ■ 9 u) — — ( du ) ■ (9 u): the second of the latter relations generates p ;/ 

while pA p combines with the contribution from 2 9p • dp, after integrating the latter 
by parts, and generates the two contributions to p'. 

From the expression for p" we see that 


dx \p"(x)\ < / dy dy' \9u(y)\ 2 \9u(y')\ / dx 


B r 


Bp X Bp 


'Br 


\x -y\\x-y’\ 


< 


(0 


2\2 


< Cp( / dy \9u(y) \ ) 

J B n 


The part with p' is more interesting: since its expression above contains inside the 
integral kernels apparently singular at x = y like \x — y\~ 1 Ap and \x — y\~ 1 9p one 
remarks that, in fact, there is no singularity because the derivatives of p vanish if 
y £ B 3p /4 (where p = 1). Hence \x — y\~ k can be bounded “dimensionally” by p~ k in 
the whole region B p /B^p/^ for all k > 0 (this remark also motivates why one should 
think p as sum of p' and p"). 

Thus replacing the (apparently) singular kernels with their dimensional bounds we get 


dx \u\\p'\ < 



dx I u 
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which can be bounded by using inequality (H) as 




where all functions depend on x (and of course on t) and then, integrating over t £ A r 
and dividing by r 2 one finds, for a suitable C > 0: 

dtdx\u\\p'\ <C(-) 1/5 G 1 r /5 Kp /5 A 1 p /5 
that is combined with I dxlullp^l < (f d^|u| 2 ) 1 / 2 (/ dx Ip"! 2 ) 1 / 2 which, in- 

V fJ -y* t/ yp yp 

tegrating over time, dividing by p 2 and using inequality (!) for / dxlp^j 2 yields: 

*J JJ r 

r~ 2 f dtdx\u\\p"\ < C(pr~ 1 ) 2 A 1 p //2 S p ). 

•J (aj r 

[3] In the context of the hint and notations for p of the preceding problem check that 
Id d® Ip'I < C(rp -1 ) 3 / dxlpl. Integrate over f the power 5/4 of this inequality, 

*J -LJ r J Dp 

rendered adimensional by dividing it by r 13 / 4 ; one gets: r ~ 13/4 LM | p '|) 5/4 < 

C(rp _1 ) 1 / 2 K p , which yields the first term of the second inequality in (4.4). Complete 
the derivation of the second of (4.4). (Idea: Note that p"(x,t) can be written, in the 
interior of B r , as p" = p -\- p with: 



1 f •£. ~ V 

P(x) = -— / -- =-p(y)u -dudy, 

J b p I— — y } 3 ~ 


9p(y) ■ (u -d)u 


p(x) — - 

4ir 


Br 


\x - y I 


dy 


(always at fixed t and not declaring explicitly the t— dependence). Hence by using 
\x — y\ > p/4, for x £ B r and y £ B p /B^ p j^, i.e. for y in the part of B p where 
9^/0) we find 
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is bounded by raising the right hand sides of the last inequalities 


to the power 5/4 and integrating over t, and finally applying inequality (H) to generate 
/ \ 5/8 

the integral ( \du\ 2 J ). 

[4]: Deduce that (4.5) holds for a pseudo-NS field ( u,p ), c.f.r. definition 1. (Idea: Let 
c p(x,t ) be a C°° function which is 1 on Q p / 2 and 0 outside Q p ; it is: 0 < p(x,t) < 1, 

\chp\ < | A(^| + |dtp < -j, if we suppose that </? has the form <p(x, t) = P> 2 ( = , ~^j) > 0 

for some <£>2 suitably fixed and smooth. Then, by applying the third of (3.1) and using 
the notations of the preceding problems, if t € A p / 2 (^o)) it is 


/ |u(x, t)| 2 dx < — / dtdx |u| 2 + / dtdx (|w| 2 + 2p) u ■ dtp < 

^Br-X{t} P Jq p Jq p 

/■ dtdx|u| 2 +| f dtdx(\u\ 2 — \u\ 2 p )u • < 9 <yJ + 2 f dtdxpu-d<p< 


dtdx\u\ 3 J +| J dtdx ^|u| 2 — \u\ 2 ju-dp H- J dtdx\p\\u\< 


<CpG p T CpJ p T p — / dfdx(|;u| 2 — |uR)u-9t/? 


We now use the following inequality, at t fixed and with the integrals over dx 


-4(L 


dx lid 


p 2 J Bp 


I U 2 - |u|2| 3 / 2 


|u| 

1 |2 1 12 
|u — uD 




~ p v ^ p 

and we also take into account inequality (P) with f — u 2 and a — 3/2 which yields 
(always at t fixed and with integrals over dx): 

/ r 1 _ _1 s/2\ 2/3 / r \ 


M - IMIJ5 




then we see that 


M 2 - Np N \dp\ < 


h\\9u\ < 


< £//2 V/2 

P 
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Integrating over t and applying (H) with exponents 3, 2, 6, respectively, on the last 
three factors of the right hand side we get 


u 


Qi 


m - m p 


< CA 1 p /2 G 1 p / 3 5 1 p /2 < C (G 2 / 3 + A P 5 P ) 


and placing this in the first of the preceding inequalities (*) we obtain the desired 
result). 


The following problems provide a guide to the proof of theorem II. Below we replace, 
unless explicitly stated the sets B r , Qr, A r introduced in definition 2, in (C) above, 
and employed in the previous problems with £>(?, A®, Q® with B ( f — {x\ \x — x_ Q \ < r}, 
A®. = {t\t 0 > t. > t 0 - r 2 }, Q° = {(x, t)\ |x — x 0 | < r, to > t > to — r 2 } = x A°. 
Likewise we shall set B ( f n = £> 3 , A® = A and we shall define new 

operators A,5,G,J,K,S by the same expressions in (3.2)%(3.8) in (C) above but 
with the just defined new meaning of the integration domains. However, to avoid 
confusion, we shall call them A°,<5°,... with a superscript 0 added. 


[5]: With the above conventions check the following inequalities 


4 u <■ n a u 
rt-n — 0/1 n+1> 


G U U 

n — n+lj 


G: 


<c« 3/2 + < 3/ V„ 3/4 ) 


(Idea: The first two are trivial consequences of the fact that the integration domains 
of the right hand sides are larger than those of the left hand sides, and the radii of 
the balls differ only by a factor 2 so that C can be chosen 2 in the first inequality and 
4 in the second. The third inequality follows from (S) with a — q — 3: 



u| 3 < C 





< C 


„3/4^0 3 / 4 




+ A 


03/2 


< 


where the integrals are over dx at t fixed; and integrating over t we estimate G® by 
applying (H) to the last integral over t .) 


[6]: Let no = n + p and Q° n — {(x, t)\\x — x Q \ < r n , to > t > t — r 2 } d = B 3 X A® 
consider the function: 


<Pn(x,t) 


exp (~(x - x Q ) 2 /A(rl + tp - t)) 
(47T (r 2 + to — t)) 3 / 2 


(x, t ) e Qn 0 


and a function Xn 0 (x,t ) = 1 on Q ° nQ _ 1 and 0 outside Q^ 0 , for instance choosing, a 

function which has the form Xn 0 (x,t) — /(r nQ x, r nQ 7 t) > 0, with <p a C°° function 
fixed once and for all. Then write (3.1) using <p = TnXn 0 and deduce the inequality 


A® + (5° 


< C 


r ~ 2 G° 2//3 + 

' n+p^n+p ' 


n+p 

E 

k=n-\-1 




0 , -2 

—I— T 

k 1 n+p 


7° 

n 


+p + 


n+p— 1 

E 

fc=n+l 


r ^ 2 Lk 


(®) 
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where L^ — r fc 2 Jqo dxdt\u\\p — p k \ with p k equals the average of p on the ball 

B ?; for each p > 0. (Idea: Consider the function <p and note that ip > (CV 3 ) -1 in 
Q 3 , which allows us to estimate from below the left hand side term in (3.1), with 
(Cr^) -1 (A® + S° n ). Moreover one checks that 


C 


C 


\<p\ < — > 1^1 < —’ n<m<n + p = n 0 , in Q^ +1 /Q m 


|dtp + A<p| < 


C 


n o 


in Q 


n 0 


and the second relation follows from dtp~\~ Ay? = 0 in the “dangerous places”, i.e. x — 
1, because p is a solution of the heat equation (backward in time). Hence the first 
term in the right hand side of (3.1) can be bounded from above by 


\u\ 2 \dtp n + Aip n \ < 


Qr 


C 


n o Jq( 


,2 C 

Ini 2 < — 

r n 0 


'Q 


|u| 3 V /3 r 5/3 < —C° 2/3 
I—I / ' n 0 - 2 ^ 

0 / 'n 0 


getting the first term in the r.h.s. of (@). 

Using here the scaling properties of the function p the second term is bounded by 


C 


n o 


\u\ \dp n \ < — 


Ini 3 + 


E § 


Qr 


n -IQ 0 


Iwl 3 < 


< 


n o 

E 

fc=n+l 


c 


fc JqI 


n-f-1 k —n+2 

n 0 


k J Ql/Ql 


k' v *fc-l 


O V—’V G9 

lnl 3 < c 22 


k=n J rl 


k 


Calling the third term (c./.r. (1.1)) Z we see that it is bounded by 


Z < 


PU ■ d_Xn 0 Pr 


o° 

^ ri( 


< 


PU ■ d_Xn+iPn 


' O 0 

^n + l 


+ 


no 

+ E . , 

fc=n+2 J Qk 

n o-l 


PHl ■ d(Xk - Xk-l)<Pr 


< 


(p- p n+1 )u- d_Xn+iPn 


' O 0 

^n +1 


+ 


+ 


E 


(p- p k )u- d(xk -Xk-l)<fr 


k=n +2 " Q°k 


+ / Ini \p\ l^(xn 0 - Xn 0 -l)pr 


Qr 


where p m denotes the average of p over B ^ (which only depends on t): the possibility 
of replacing p by p— p in the integrals is simply due to the fact that the 0 divergence of 
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u allows us to add to p any constant because, by integration by parts, it will contribute 
0 to the value of the integral. 

From the last inequality it follows 



+ J° r~ 2 
' «0 n 0 


fc=n+1 


n 0 — 1 

E 

fc=n +1 


J L k 

r k 


4- J° r ~ 2 
' n o n 0 


then sum the above estimates.) 


[ 7 ] If x Q is the center of Q the function Xn 0 P can be regarded, if no < — 1 , as defined 
on the whole R 3 and zero outside the torus Q. Then if A is the Laplace operator on 
the whole R 3 note that the expression of p in terms of u ( c.f.r . (a) of (3.1)) implies 
that in Q® : 

Xn 0 P = A -1 A Xn 0 P = A -1 (pAxno + 2 (dxn 0 ) ■ (dp) ~ Xn 0 dd ■ (uuU 
Show that this expression can be rewritten, for n < no, as 

Xn 0 p = - A ~ i (Xn 0 dd(uu)) + [A -1 (pAXn 0 ) + 2(dA~ 1 )((dXn 0 )p)] = 


= [-(ddA 1 )(xn 0 uu)] + [2(0A L )(dxn 0 uu ) - A 1 ((ddXn 0 )«u)] + 


-l 


-D 


+ [A 1 (pAxn 0 ) + 2(dA 1 )((dxn 0 )p)] d = PI + P2 + P3 + P4 
with pi = -(ddA~ 1 )(xn 0 '^n+iuu) and p 2 = -(ddA~ l )(xn 0 (l ~ $n+i)uu ) where 

dfc is the characteristic function of B°^ and P 3 ,P 4 are the last two terms in square 
brackets. (Idea: Use, for x,t £ Q® , Poisson formula 


Xn 0 (X, t)p(x, t ) = 


-1 

47T 


A ((xn 0 p)(y,t)) 


dy — 


’ B° 

n 

-l 

4n j 

and suitably integrate by parts). 


\x-y\ 

pAxn 0 + 2 d_Xn 0 ■ dp - Xn 0 dd ■ ( uu ) 


si 


|x - y | 


■dy 


[8] In the context of the previous problem check that the formulae derived there can 
be written more explicitly as 


Pi=~(ddA x ) • (xn 0 Vn+iuu), 


P 2 = - — 
4tt 


f B°/BO n+1 


El) ' 


u 


I x-y 




’ B 


x — y 1 

■(9xn 0 )uu + — 


0 \x-y\ 




|x - y| 


(§dxn 0 )uu 


PA 


177 


1 2 

]-rP(y)Axn 0 + — 

\x — y\ — 47r 


x — y 

p(y ) I-ti 

— |x — y\ a 
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where n < uq and the integrals are over y at t fixed, and the functions in the left hand 
side are evaluated in x,t. 

[ 9 ]: Consider the quantity L n , introduced in [6], 


def _ 2 


L n = r 



u\ |p 


Pn(i9)\d^cM 


and show that, setting uq — n + p, p > 0, it is 


L n <C 

V r n +1 > 

7/5 , 

\ A 0 V+0 1 / 5 7>^0 4 /5 , 

1 ^n + 1 ^n+1 A n 0 + 


L \ 

i n 0 / 

n 0 -l 


+ G n+1 

1 r 3 r 0 V3 r~ 3 A° 

' r n+l'^ r n+l / J r k 





fc=n+l 



r n +1 


n 0 


5/3 


V+0 2 / 3 , 

< - T n+1 '- T n o ^ 


(Idea: Refer to [8] to bound L n by: y~\_ T r n 2 J^ 0 |u| |pi — p ™| where p™ is the average 

of pi over £>0; and estimate separately the four terms. For the first it is not necessary 
to subtract the average and the difference |pi — p™ | can be divided into the sum of 
the absolute values each of which contributes equally to the final estimate which is 
obtained via the (CZ), and the (H) 


\pi - Pi\\u\ < 2 


bi| 3/2 


’K+i 


’ B ° n+ i 


2/3 


u 


'B° 

n+1 


1/3 


< c 


u 


' B° 

n + 1 


and the contribution of pi at L n is bounded, therefore, by + note that this 


would not be true with p instead of pi because in the right hand side there would be 
J" \p\ 3/2 rather than J 0 |p| 3 / 2 , because the (CZ) is a “nonlocal” inequality. The 


n + 1 


term with p 2 is bounded as 


J A°„° JB° 

< r n 


\P 2 ~ P2 llltl < 


A 


00 


|u| r n max \dp 2 | < 




u 


3 \ 1/3 


2/3 10/3 

r ' rr* ' 

n. ' 71 . 


Qn U 
no~ 1 


max \dp2 1 < 
Qn 


< r 5 G° 1/3 

— 1 n^n 


E 

m=n-\-1 


max 

te A° 


u 




1 /, A° 


= r b G' 

1 n 


E + 

rpCj 


m.=n-\-l 
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Analogously the term with ps is bounded by using \dp 3 \ < Crj/ 4 f 0 |u| 2 which is 

B n 0 

majorized by CY^ 3 ( f B0 |u| 3 ) 2 / 3 obtaining 


\u\\p 3 - Pl I < -^-D, 3 / [( / |u| 3 ) 2 / 3 ry 


lull < 


< £ r 3 r -3 / ( / 

- r 2 n n 0 / W 
J A° 


TC d B o 


(/ N 3 ) 2/3 ( / N 3 ) 1/3 < 


< ^( —) 3 dfr; / 3 G»„ 2 / 3 G» 1/3 = C(-^) 5 / 3 G“ 0 2 / 3 GS 1/3 


2 v 

/y» Zi ry* 

' n r n 0 


Finally the term with p^ is bounded (taking into account that the derivatives A Xn , d_Xn 
vanish where the kernels become bigger than what suggested by their dimension) by 
noting that 


/ |P4 - P? | H < Cr n / |u| max 1dp 4 1 < Cr n ( / |u| ) ( / —— ) 

J B o JbO bo \J b o J\J B O o < 0 J 

Denoting with the operator iC® 0 without the factor r nQ 13 ^ 4 which makes it di¬ 

mensionless, and introducing, similarly, A®, (5® we obtain the following chain of in¬ 
equalities, using repeatedly (H) 

:/' ^ o, U(f (f A YT(f (( 4 T-- 


r n A 0 ' JB° r n 0 

n n 0 


A 0 k Jb° 

^n n 


5\ 1/5 




' A° x Jb° 

n n 


^< 4/5 ( 

' n n 0 v J B° 


l 3 j 1/5 d /5 < 


, C r n 12 / 5 - 0 4/5 - 0 1/5-0 1/5 / r n 

< ~~o 4 r n K n 0 A n G n < G [ - 

ry* r-t ry* U \ ry 


nr* Zj rv*^t 

' n ' n o 


40 1/5^0 1/5 t>-0 4 /5 

^ n ^ tiq 


Finally use the inequalities of [5] and combine the estimates above on the terms pj , j = 

1, ••, 4.) 

[10] Let T n = (A® + h 3 ); combine inequalities of [6] and [9], and [5] to deduce 

n 0 — 1 n Q -1 

T n < 2 2n (^2~ 2no £+ ^ 2 _2fc T 3/2 + 2 _2n °£ + 2 _7n ° /5 e ^ 2 _3fe / 5 T^ /2 + 

k=n-\- 1 fc=n+l 

rig — 1 ng —1 no — 1 np-1 

+ £ 2 _5n °/3 2 _fc/3 T fc 1/2 + 2 _2fc T 3/2 + ^ 2 fc T fc 1/2 2~ 3p r p ^] 

fc=no+2 fc=n+l fc=n-|-l p=fc 


e = C'max(G° 2/3 ,iF° 4/5 , J° no ) 
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and show that, by induction, if e is small enough then r n 2 T n < e 2 ^r n 2 . 

[ 11 ] : If G(ro ) + J(ro) + K(tq) < e s with e s small enough, then given (x',t') G 
Qvq/4 Olio ) ^0) j show that if one calls G®, J®, X®, A®, <5® the operators associated with 
Qr(x' ,t') then 

1 2/3 

limsup — A® < C £s 9 

rvt Z Z 

n^oo 1 n 1 0 

for a suitable constant C. (Idea: Note that Q® o / 4 (a/,tl) C QroGkcu^o) hence 
G® o / 4 , J° o y 4 , . . . are bounded by a constant, (< 4 2 ), times G(r o), J(ro). . . respectively. 
Then apply the result of [10]). 

[ 12 ] : Check that the result of [11] implies theorem II. (Idea: Indeed 



\ A n > \ I I u(x,t')\ 2 dx- 

rpA rpO 

1 n 1 n 


4-7T , . . 

—— \ u (x , t) 

> — OO o I — k— ' / 


where 13® is the ball centered at x', for almost all the points (x/ ,t r ) G Q® q / 4 ; hence 
\u(x',t')\ is bounded in Q® o / 4 and one can apply proposition 2). 

[13]: Let / be a function with zero average over 13,9. Since /(x) = f(y) + f 1 ds <9/(y + 
(x — y)s) • (x — y ) for each y G 13®, averaging this identity over y one gets 


/(t) 



dsdf(y + (x 


y)s) • (x 


y) 


Assuming a = 1 prove (P). (Idea: Change variables as y —► z — y + (x — y)s so that 
for a integer 



I/(t)I 


a dx 
|13,0| 




dz_ 

Mi a 


Sf(z) ■ (z 

— s) d 


x) 


a 


where the integration domain of z_ depends from x and s, and it is contained in the 
ball with radius 2(1 — s)r around x. The integral can then be bounded by 


dz, ds i 


dz„, ds c 


113,® | 1 — si 113® | 1 — 


■(2rr\df(z 1 )\...\df(zJ\ 


dx 

MI 


where x varies in a domain with |x — zA < 2(1 — Si)r for each i. Hence the integral 


dx 


over is bounded by 8(1 — s^) 3 for each i. Performing a geometric average of such 
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bounds (over a terms) 


B 


0 l/( - )r ra 


< 2 a + 3 r a 


CX. n 

n/ 


d - idS ' l£tf(ii)l(l-Si) 3/ “< 


2—1 
OL 


Oi+2>n 


< 2 r 


(J B „ mu)l w\) i 


1 


ds 


(1 — s ) 3 3 / c 


getting (P) for a = 1 and an explicit estimate of the constant Cf: this also gives a 
heuristic motivation for (P) with a < ■|.) 

[14] : Differentiate twice with respect to a 1 and check the convexity of a 1 —> 
II/Hq, = ( J\f(x)\ a dx/lB®]) 1 / a ■ Use this to get (P) for each 1 < a < ao if it is 
valid for a = ao. (Idea: Since (P) can be written ||/|| a < Ca (f \df\dx/r 2 ) then if 

a -1 = da^ 1 + (1 — d)(ao + 1) _1 with ao integer it follows that C a can be taken 
Col — dCoLQ + (1 — 1 ?) Cct Q -\-l). 

[15] : Consider a sequence u ^ of solutions of the Leray regularized equations which 

converges weakly (i.e. for each Fourier component) to a Leray solution. By construc¬ 
tion the verify the a priori bounds in (0.2) and (hence) (2.5). Deduce that u 

verifies the (3.1). (Idea: Only (c) has to be proved. Note that if —> u° weakly, then 
the left hand side of ( 1 . 1 ) is semi continuous hence the value computed with uP is not 
larger than the limit of the right hand side in (1.1). On the other hand the right hand 
side of (1.1) is continuous in the limit A —> oo. Indeed given N > 0 weak convergence 
implies 


lim 

A—>oo 






< lim 

A—>oo 


< lim 

A—>oo 


( E 

0<|fc|<AT 

( E 

0<Ifcl<AT' 


° (t)| 2 + /’ ° d t%h x k m - 2 » ( t)| 2 ) 


< 


I fc I > AT ' 


'T 0 


'T 0 




dt / L§(« A - « 0 )| 2 ) 


’n 


'To 


= lim —— 

A—>oo 


dt 


[§(h a -«°)l 2 < 


2Eqis 


-l 


N 2 


using the a priori bound in ( 0 . 2 ) (with zero force) and componentwise convergence 
of the Fourier transform 7 ^ (t) of u(t ) to the Fourier transform 7 ® (t ) of uP. Hence 

0 showing the convergence of the first two terms of the right 
hand side of (1.1) to the corresponding terms of (c) in (3.1). 
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Apply, next, the inequality (S), (2.2), with q — 3, a — |, a ~ \i an< ^ a § a i n by the 
a priori bounds in (0.2) we get 



r T <o 

< C(2S 0 ^- 1 ) 3/4 (2 V /^) J dt\\u x -u°\\ 2 2 ^^0 


showing continuity of the third term in the second member of (3.1). Finally, and anal¬ 
ogously, if we recall that p x — —A -1 ^ didj(u x u x ) and if we apply the inequalities 
(CZ) and (H), we get 


dt dx\p x u x —p°u° | < J J \p x — p°\ |u A | + J j |p°| | u x — u°| < 

( f J \p x - P° |3/ 2 ) 2/3 (/ / I U X I 3 ) 1/3 + ( / / \p° | 3/2 ) 2/3 (/ / h X - U° I 3 ) 1/3 


where the last integral tends to zero by the previous relation while the first, via (CZ), 

will be such that 0 \p x — p° | 3 / 2 < ( J J |u A — u°| 3 J A _ oo > 0 proving the 

continuity of the fourth term in the right hand side of (c) in (3.1). Hence the right 
hand side is continuous in the considered limit). 

[16] : (covering theorem , (Vitali)) Let S be an arbitrary set inside a sphere of R n . 

Consider a covering of S with little open balls with the Vitali property: i. e. such that 
every point of S is contained in a family of open balls of the covering whose radii 
have a zero greatest lower bound. Given p > 0 show that if A > 1 is large enough 
it is possible to find a denumerable family F\, F 2 , ■ ■ ■ of pairwise disjoint balls of the 
covering with diameter < p such that U 7 ;Aiq D S where AF) denotes the ball with the 
same center of Fi and radius A times longer. Furthermore A can be chosen independent 
of S, see also problem [17]. (Idea: Let T be the covering and let a — maxjr diam(F). 
Define a^ = a2~ k and let J~\ be a maximal family of pairwise disjoint ball of T with 
radii > a2“ J and < a. Likewise let IF 2 be a maximal set of balls of T with radii 
between a2 —2 and a2 -1 pairwise disjoint between themselves and with the ones of 
the family T\. Inductively we define 3 ~\, . . ., Tq,, .... It is now important to note that 
if x 0 U kFk it must be: distance(x, F^) < Xa2~ k ' for some k, if A is large enough. 
If indeed 5 is the radius of a ball S§ containing x and if a2~ k ° < 5 < a2 _/ '° + 1 then 
the point of S§ farthest away from x is at most at distance < 25 < 4a2 —fc °; and if, 
therefore, it was d(x,Fk 0 ) ^ 4a2 _fc ° we would find that the set Fk 0 could be made 
larger by adding to it S§, against the maximality supposed for Fk 0 - Note that A = 5 
is a possible choice.) 

[17] : Show that if the balls in problem [16] are replaced by the parabolic cylinders 
which are Cartesian products of a radius r ball in the first k coordinates and one 
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of radius r a , with a > 1 in the n — k remaining ones, then the result still holds if 
one replaces 5Fj with A Fi where A is a suitable homothety factor (with respect to 
the center of Fi). Show that if a = 1,2 then A — 5 is enough (and, in general, 
A = (4 2 + 2 2 ( 1 +“)/ a ) 1 / 2 is enough). 

[18]: Check that the Hausdorff dimension of the Cantor set C is log 3 2, c.f.r. (5.3). 
(Idea: It remains to see, given the equation in footnote 9 , that if a < ao then /a a (C) = 
oo. If 5 = 3~ n the covering C n of C with the n—th generation intervals is “the best” 
among those with sets of diameter < 3~ n because another covering could be refined 
by deleting from each if its intervals the points that are out of the n —th generation 
intervals. Furthermore the inequality 1 < 23 -a for a < log 3 2 shows that it will not be 
convenient to further subdivide the intervals of C n for the purpose of diminishing the 
sum \Fi | a . Hence for S — 3~ n the minimum value of the sum is 2 n 3~ na n _ >00 > oo.) 
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